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ABSTRACT 

We present an analysis of the ionic composition of iron for two interplanetary 
coronal mass ejections observed in May 21-23 2007 by the ACE and STEREO 
spacecraft in the context of the magnetic structure of the ejecta flux rope, sheath 
region, and surrounding solar wind flow. This analysis is made possible due to 
recent advances in multispacecraft data interpolation, reconstruction, and visu- 
alization as well as results from recent modeling of ionic charge states in MHD 
simulations of magnetic breakout and flux cancellation CME initiation. We use 



Of 



these advances to interpret specific features of the ICME plasma composition 
resulting from the magnetic topology and evolution of the CME. We find that 
in both the data and our MHD simulations, the flux ropes centers are relatively 
cool, while charge state enhancements surround and trail the flux ropes. The 
magnetic orientation of the ICMEs are suggestive of magnetic breakout-like re- 
connection during the eruption process, which could explain the spatial location 
of the observed iron enhancements just outside the traditional flux rope magnetic 
signatures and between the two ICMEs. Detailed comparisons between the simu- 
lations and data were more complicated, but a sharp increase in high iron charge 
states in the ACE and STEREO-A data during the second flux rope corresponds 
well to similar features in the flux cancellation results. We discuss the prospects 
of this integrated in-situ data analysis and modeling approach to advancing our 
understanding of the unified CME-to-ICME evolution. 
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Subject headings: Sun: solar-terrestrial relations — Sun: coronal mass ejections 
(CMEs) — Sun: flares — Sun: heliosphere — Sun: solar wind 



Introduction 



Coronal mass ejections (CMEs) are periods of explosive magnetic energy release and 
coronal field reconfiguration that blow out huge portions of the quasi-stable solar atmo- 
sphere into interplanetary space (jTouseyl 1 19731 : iGoslind Il997l) . When these solar eruptions 
are Earth-directed, they generally arrive in 1-5 days (ICane &: Richardson! 120031) and thei r 



impact can cause geoeffective space weather responses (IGosling 



1993 



Zhang et all 120041 1 . 



CMEs are observed soon after their eruption by remote sensing instruments that capture 
the global structure of these events. Their morphological structure in coronagraph obser- 
vations typically consists of a bright bubble of plasma, sometime s with a darker cavity and 



bright core (e.g., Illling fc Hundhausenlll986t iHoward et al.l 119971 1 . Once these CMEs prop- 



agate into the heliosphere they are commonly referred to as interplanetary CMEs (ICMEs) 
and are observed in-situ by spacecraft that sample local measurements of the magnetic 
field, plasma, and composition. ICME identification is based on a variety of parameters 
inc luding enhancements in magne tic field, composition and temperature depressions (e.g., 



sec 



Zurbuchen fc Richardson! 120061 . and references therein) 



Though CMEs and ICMEs are the same phenomena, the remote sensing and in-situ 
observations can often be difficult to directly relate to one another, due to both the na- 
ture of the observations as well as the evolution that takes place during the transit from the 
Sun to the spacecraft. In-situ solar wind composition data provides an important connection 
between the ICME observations and associated observations at the Sun, in particular provid- 
ing a measure of physical properties in the corona that can be compared to EUV and X-ray 
emission and white light coronagraph observations. Charge state composition is directly 
related to the electron temperature at the source region , with hotter source material result- 
ing in higher in-situ charge states (jBurgi &: Geissl Il986h . During the solar wind expansion, 
these char ge states freeze-in once t he expansion timescale becomes larger than the collision 
timescale (IHundhausen et al.lll968l ). ICMEs often contain unusual charge state distributions. 
These charge states can be enhanced relative to the ambient solar wind due to heating at 
the source region, lower than the ambient solar wind due to incompletely heated cold promi- 



nence material, or a combination of both low and high charge states (IBame et al 



Gloeckler et al. 



2011 



1999 



Gilbert et al. 



Richardson &: Canell2004t iLepri fc Zurbuchen! |2010| ; 



1979; 



Gruesbeck et al. 



20121 ) . Because of this complexity, it is important to examine the full 
charge state distribution. Charge state enhancements are weakly correlated with associ- 
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ated flare mag nitude (jReinardll2005l ) and tend to be strongest near the center of the ICME 
( Reinard1l2008 ). In this paper we examine the in-situ iron charge state data from the ACE and 
STEREO spacecraft and compare them with theoretical predictions derived fro m MHD simu- 



l ation s of two different CME initiation models using the technique developed by lLynch et al. 

(bonh . 

Given the localized nature of ICME observations (ICMEs are typically observed by a sin- 
gle spacecraft that samples a narrow trajectory through the larger CME structure) drawing 
conclusions about the global structure of ICMEs can be complicated as it is not always clear 
which part of the ICME is being sampled. The subset of ICMEs with the most structured field 
and plasma signatures are known as magnetic clouds and are defined as having an enhanced 
field strength, relatively s mooth magnetic field rotations, and a lower than expected tempera- 
ture ( IBurlaga et al.lll98ll ). While the co mmonly accepted occurrence rate for ma gnetic clouds 
is between 30-50% of all ICME events (lGoslind[l99ol ; iRichardson k Cane! 120041 ) . there is in- 



(Jian et al. 


2006; 


Kilpua et al. 


2011) 



criteria to be considered magnetic clouds (IBurlaga et al.lll98ll ) the overall ICME structure 
can be approxi mated by modeling t he in-situ ICME magnetic fields as force-free, cylindrical 
flux ropes (e.g.. lLepping et al.lll990f ). These techniques allow researchers to interpret the cav- 
i ty portion of the CME coronagraph morphology with the in-situ magnetic flux rope stru cture 



(Illine k Hundhausen 


1986; 


Wood et al. 


1999: 


Webb et al. 


2009; 


Voulidas et al. 


2012 


) and 
vever, 


indeed observations have confirmed this association (e.g. 


DeForest et al. 


2011 


). Hov 



analysis of the current density structure, measures of the interplanetary pressure-balance, 
and multispacecraft observations have shown that even magnetic clouds are not precisely 
force- free, cylindrical ob jects. Therefore non-force free and consequently non-cylindrical 



fits w ere developed (e.g.. iFarrugia et al.lll993l . Il999t iMulligan k Russell! \2 00 It iHidalgo et al. 



20021 ) that allowed more complexity in the local ICME field modeling, while still retaining a 



big picture view of the global structure. Though these advances have allowed a greater num- 
ber of ICMEs to be modeled, all of these analytic flux rope modeling techniques are limited in 
that they are only relevant for the flux rope-like portion of the ICME; non-flux rope ICMEs or 
cuts through an ICME that do not cleanly intersect the flux rope cannot be included because 
they lack the proper magnetic structure needed for the fits. In theory, t he pr e ssure-balance 



Grad -Shafranov reconstruction techniques (e.g., iHu k Sonnerupl l200ll . |2002| ; llsavnin et al. 



20 111 ) avoid the a-priori flux rope criteria, as well as the problematic identification of the 



flux ro pe boundaries, but in practice are most su ccessful wh en reconstructing flux rope struc- 
tures flRilev et al.lbooi kl-Haddad et al.lhoilh . Recently. [Mulligan et all feoi2h addressed 
this issue by developing a technique to incorporate non-MC (magnetic cloud) material into 
MC flux rope modeling. 
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Multispacecraft observations of ICMEs have made and continue to make a valuable con- 
tribution to our understanding of the larger heliospheric structure, placing important con- 
straints on flux rope modeli ng; and the spatial extent of ICME ejecta over which the coherent 



flux rope signatures persi s t (IBurlaga et al. 



19811 : lHammond et al.lll99a : Mulligan et al.lll999 



Mulligan fc Russell! U 00 ll; iRilev et al.ll2003h . The Solar-Terrestrial Observatory (STEREO) 
mission (IKaiser et al.ll2008l ) was launched in October 2006 and consists of twin spacecraft at 
1 AU traveling ahead of and behind the Earth at a rate of approximately 22.5 degrees per 
year. The STEREO observing geometry provides a means of relating multispacecraft in-situ 
observations to each other and to the continuous remote coronal and helios pheric imaging 



start ing at the Sun and extending t o 1 AU. The data from the SECCHI suite (Howarc 



20081 ) combined with the IMPACT flLuhmann et al.ll2008f ) and PLASTIC flGalvin et al. 



et a 



20081 ) 



in-situ field and plasma measurements have ushered in a new era of direct CME-ICME obser- 



vations, explicitly linking various morphologica 



ties (e.g., 


Harrison et al. 


2008, 


2010; 


Davis et al. 


2009 


Davies et al. 


2009; 


M 


ostl et al. 


2009a: 


Roulliard et al. 


2009 




Webb et al. 


I2009 




Lvnch et al. 


2010: 


DeForest et al. 


2011 


). In addi- 



features to their associated in-situ proper- 



tion, the two STEREO spacecraft, along with the ACE and WIND spacecraft, offer three or 
four tracks through a give n ICME rather tha n just one, allowing more of the ICME spatial 
structure to be sampled (jKilpua et al.ll201ll ). However, often even three tracks of obser- 
vations through an ICME cannot fully constrain the ICME global structure and further 
heliospheric modeling is necessa ry to interpret and re late the multispacecraft data back to 
remote observations at the Sun. iMulligan et al.l ( 120121 ) describe a sophisticated method that 
provides an interpolated spatial mapping of the data between the spacecraft tracks. This 
method provides an estimate of the global ICME structure and simplifies the comparison 
with modeled and observed CME structures near the Sun. 

The structure of the paper is as follows. In Section 2 we describe our methodology. First, 
in section 2.1, we briefly review previous work on the source region and the multispacecraft 
i n-situ data of the 2007 May 21-23 ICMEs. In section 2.2 we describe the IMulligan et aL 
(120121 ) reconstruction and interpolation techniques for transforming the in-situ timeseries 
data into spatial maps and present two-dimensional composition data in the context of the 
ICME flux rope st ructures and their surrounding regions. In section 2.3 we describe the 
Lynch et al.l ( 1201 ll ) application of ionic charge state calculations to numerical MHD simula- 
tions to derive theoretical spatial distributions of iron charge states for two idealized CME 
initiation scenarios. In section 3 we compare the observed and simulated ionic composition 
structure; first, on the global scale of the two flux rope ICMEs and their spatial structure, 
and second, by constructing a set of synthetic spacecraft trajectories through the MHD sim- 
ulation data for direct comparison with the multispacecraft measurements of the detailed 
iron charge state distributions (Fe +6 to Fe +20 ). In section 4 we discuss how the structure 
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and complexity of the interplanetary composition data is reflected in the complexity of the 
May 19 and 20 CME source region's magnetic topology and resulting eruption scenario. We 
conclude, in section 5, with a summary of our results and prospects for future analyses. 



2. Methodology, Data Analysis, and Numerical Modeling 
2.1. Event Description and Multispacecraft Observations 



The 2007 May 21-23 ICMEs occurred when STEREO-A and STEREO-B were approx- 
imately 9 degrees apart with a separation from ACE of 5.9 and 3.1 degrees, respectively. 
These events both originated in AR10956, which was located nearly exactly at disk center. 
The first event was associated with a B9.5 flare that peaked at 13:02 UT on May 19, while the 
second event was associated w ith a B6. 7 flare that peaked at 05:56UT on May 20. The first 
of these events is well studied. iLi et al.l (120081 ) examined the sourc e region topology and the 
photospheric magnetic field evolution leading up to the eruption . iBone et al.l (120091 ) exam- 
ined the pre-emption filament formation and lLiewer et al.l (120091 ) an alyzed the stereoscopi c 
3-dimensional traj e ctory durin g the fi lamen t lift-off and eruption. IVeronig et al.l ( 120081 ) 



Gopalswamy et a 



wave signatures. 



Kerdraon et al 



(2009) and lAttrilll (120 lOh characterized the EUV-dimming and coronal 



(120101 ) examined the CME-ICME connection in terms of 
the local acceleration of energetic electrons during the eruption and their interplanetary ra- 
dio burst signatures. The left column of Figure [1] shows the 17lA EUV dimming associated 
with the May 19 eruption. In the top panel is the base image, taken at 12:22UT, in which 
the active region is the bright region seen in the c enter. In the bottom panel i s the difference 
image (13:12UT-12:22UT, technique described in lReinard fc Biesecker 1120081 ) with the dark 
dimming regio n indicating the erupting materia l from the northwest section of the active 
region (see also lLi et al.ll2008l ; iKilpua et al.ll2009l ). The right column in Figure [T] shows sim- 
ilar images for the May 20 eruption (04:52UT base image, top panel, and 06:02UT-04:52UT 
difference image, botto m panel). In this ca se, the eruption originated from the southeast of 
the same active region. iMierla et al.l (120081 ) utilized STEREO/COR1 height-time diagrams 
to analyze the propagation direction of the May 20 eruption. 



Kilpua et al.l (120091 ) discussed the CME-ICME connection during this period associating 



the CDAW catalog entry of a fast CME occurring at 13:24 UT (speed of 958 km/s) with 
the well-structured magnetic cloud observed by STEREO, ACE, and WIND on May 22. 
A second, slower C ME was observed at 13:48UT (speed of 294 km/s) and associated by 



Kilpua et al. 



Kilpua et al. 



(2009) with a second magnetic cloud seen by STEREO-A on May 23. Following 
J2009h we denote these events as MCI and MC2. Figure [2] shows the STEREO- 



B, ACE, and STEREO-A in-situ data (from left to right) for the 2007 May 21-23 period 
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(DOY 141-143), which includes ICME signatures resulting from both eruptions. From top 
to bottom, we plot the bulk radial proton velocity v p , proton number d ensity n p , pr oton 
temperature T p (black line) and the expected temperature T exp (green line; lLopezlll987l ). the 
distribution of iron charge states (Fe +6 to Fe +20 , denoted as QFe), and the interplanetary- 
magnetic field in RTN coordinates: -Br (red), B<? (green), (blue), and B magnitude 
(black). In each panel the MCI interval is indicated by the yellow vertical lines. Spacecraft 
data that clearly transverse the magnetic flux rope portion of the ICME are shaded. The 
MC2 interval is likewise indicated in light blue lines and shading. 



Both lLiu et all (120081 ) and lKilpua et al.l (120091 ) performed Grad-Shafranov fits for MCI 
obtaining moderate ly inclined flux rope orientations with reasonably circular cross-sections. 
Mostl et al.l ( l2009bl ) obtained similar re sults applying a mu ltispacecraft Grad-Shafranov fit 
to STEREO-B and WIND data, while iMostl et all (bo09ch fit MC2 using Grad-Shafranov 
modeling at STEREO-A, which was then compared and optimized using WIND observations 
to yield a distorted, slightly flattened, ellipse. 



Rakowski et al.l (120111 ) examined the ionic charge states during MCI and modeled the 



evolution of the flux rope ejecta with an analytic spheromak magnetic field configuration. 
They calculated a predicted ionic charge state distribution resulting from internal ICME 
heating supplied via the dissipation of magnetic energy associated with the ICME magnetic 
structure relaxing toward a force-free equilibrium. This approach was able to produce Fe +8 
to Fe +13 as well as other moderate charges states in C, O, and Si, but not t he observed 
shapes of the full charge state distributions. We describe the iRakowski et al.l ( 120071 . l201ll ) 
procedure for the calculation of ionic charge states and our application of it to our MHD 
simulations of CME initiation and flux rope formation in Section 2.3. 



2.2. Spatial Mapping and In-Situ Data Interpolation 



Mulligan et al.l (120121 ) have recently develop ed important extensions to multispacecraft 



data analysis techniques. First, by extending the lMulligan &: Russelll ( 120011 ) non-cylindrical 
flux rope model to include the bulk plasma velocity vector into the model inversion and 
fitting, and second, by using the velocity time series in multispacecraft data to create inter- 
polated spatial maps of the solar wind and ICME plasma structures. By incorporating the 
velocity flow deflections around the ICME obstacle into the calculation of the orientation 
and size of the flux rope portion of the ICME, we can utilize multispacecraft data clearly 
associated with the envelope of the ICME but not the flux rope itself in both the ICME 
modeling and the interpretation of the surrounding heliospheric solar wind structure. The 
spatial mapping procedure consists of first integrating the components of the bulk velocity 
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time series in each of the multispacecraft data tracks to obtain the relative spatial positions 
of any measured or calculated in-situ quantities of interest and then empl oying Delauney 



triangulation to interpolate between each of these spacecraft tracks. The iMulligan et al. 



( 120 12l ) spatial mapping procedure i s derived from interpolation techniques that have been 
used in the analysis of Cluster data (jChanteur fc Motteall993l ; IChanteurlll998l ; ISibsorull98lh 
and rep resents a sub s tantia l imp rovement over the sim pl ified version of t he spa tial mapping 
used in iLynch et al.l ( 120031 ) and iReinard et al.l (120101 ) . IMulligan et al.l (120121 ) constructed 
spatial contour maps from the STEREO, ACE, and WIND time series data during the 2007 
May 21-24 period to analyze the ICME and surrounding solar wind structure of MCI and 
MC2 in the bulk plasma velocity, proton density and temperature and He ++ /H + ratio (see 
their Figure 9). 

In Figure E] we plot the spatial mapping of the ACE/SWICS and STEREO/PLASTIC 
ionic composition measurements of the average iron charge state (QF e )- In this plot the time 
series data is converted to spatial coordinates with the origin of the spatial map located 
at the Earth's position in the y dire ction and centered between the two flux ropes in the x 
direction (as in IMulligan et al.ll2012l ). The axes are given in AU and the individual spacecraft 



tracks are shown in grey, with STEREO-B (STB) tracing the bottom edge of the spatial map, 
ACE in the center, and STEREO-A (STA) tracing the top edge. Each spacecraft's position 
in the plot increases from right to left in time and the Sun is to the left. The black arrow 
vector field regions are the GSE (Geocentric Solar Ecliptic) B x and B y components of the 
interpolated magnetic field showing the large-scale rotation associated with the cross-sections 
of the MCI and MC2 flux rope structures. Specifically, the vector arrows drawn along the 
spacecraft track are the observed data in the yellow (blue) intervals shown in Figure [2] for 



MCI (MC2) and the arrows between spacecraft tracks are determined by the IMulligan et al. 



( 120121 ) interpolation method. From the spatial mapping of the flux rope structure it is clear 
that STB passes near the center of MCI and is either at or just outside the boundary of 
the MC2 flux rope; ACE passes comfortably through both flux ropes, but has a slightly 
higher impact intersection at MCI with respect to the "center" of the cross-sectional field 
rotation; STA just misses th e edge of MCI and has a relatively high-impact angle intersection 
with MC2. As discussed by IMulligan et al.l ( 120121 ). the flux rope structures inferred by the 
spatial mapping interpolation are completely consiste nt with various multispacecraft flux 



rope modeling res ults obtained by other researchers ( jLiu et al.l l2008t iKilpua et al.l 12009 
Mostl et al.ll2009 



Examining Figure E] we see low values of (QFe) m the front (right) portion of MCI and 
at the flank that crosses STA. There is an increase in (Q-F e ) in the back half of MCI along 
the ACE trace through the data. Interestingly, STB, which passes most closely through the 
center of the flux rope, does not detect any (Qf c ) enhancements at all. For MC2 the center 
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of the flux rope has (Qf c ) values that are significantly lower than the rest of the flux rope. 
The highest values of (Qp e ) in MC2 occur at the STA flank and following the event. Also of 
interest is the (Q-Fe) enhancement that seems to link the two events, extending from the back 
half of MCI, filling the space in between the two events, and blending into the enhancement 
at the STA flank of MC2. 

We have used (Qfc) in the spatial mapping because i t is a good proxy for i ron ch arge state 



enhancements and indicative of ICME material (e.g., iLepri fc Zurbuchenl 120041 ). however 
there are limitations to using just the distribution average alone. By comparing the timeseries 
of the QFe distributions shown in Figure 2 with the Figure 3 spatial map, we see that there 
are two different signatures that give rise to similar enhancements in (Qf c )- For example, the 
ACE MCI flux rope and STA data preceding the MC2 flux rope have a (QFe) enhancement 
due to a general broadening of the entire distribution to moderately higher charge states, 
whereas the STA MC2 (QFe) enhancement arises because of a bimodal QFe signature that 
includes a significant peak at the high iron charge states of +15/+16. 



2.3. Inferring ICME Charge States from MHD Modeling 



Lynch et al.l (120111 ) demonstrated a post-processing procedure in which the ionic charge 
state composition of heavy ions is calculated from MHD simulations of CME initiation using 
the single-fluid (proton) velocity, temperature and density evolution. Starting with a snap- 
shot late in the modeling run, a set of individual fluid parcels are integrated backwards in time 
using the MHD velocity field to derive the temperature and density history of these parcels. 
The MHD temperature and density profiles are then used as input s into the Blast Prop 



agation in a Highly Emitting Environment co de [BLASPHEMER] lLaming fc Grunl 12002 



Laming Sz Hwang! 120031 ; iRakowski et al.l 120071 ) which solves the continuity equations for 
heavy ion charge states of interest, 



dn, 



Ot 



i - {c'r + <T + cf ) n q + (c^ + C* J n q+1 ] 



where n e is electron density and n q is the density of a given charge state. Here, for a 
charge state q, the right hand side represents the sources and sinks of this charge state from 
electron impact ionization (C lon ), and radiative (C rr ) and dielectronic (C dr ) recombination. 
This approach makes the following assumptions: (1) the heavy ions act as passive tracers 



of the MHD flow, (2), the plasma remains electrically neutral n e 



n,. 



(where p refers 



Lynch et all (1201 lh 



to protons), and (3), the electron temperature is given by T e = T p 
applied this technique to two different axisymmetric MHD simulations run with diff erent 
codes; a magnetic breakout CME run with the ARC7 code (lAllred fc MacNeicd 120081 ) and 



-9 - 



a flux cancellation CME run with th e Magnetohydrodynamics-on-A-Sphere (MAS) code 
(ILinker et al.ll2001t iLionello et al.ll2009l ). The comparison of the derived charge states of the 
ARC7 magnetic breakout simulation and the MAS flux cancellation simulation provided 
strong evidence that the eruptive flare heating in the low corona could produce enhancements 
in commonly measured high charge state ratios (e.g., 7+ /0 6+ , Fe- 16+ /Fetotai ) within the 



CME flux rope that would be carried out into the heliosphere. As mentioned in lLynch et al. 



( 120111 ). these idealized simulations were not designed for the ionic charge state analysis we 
performed or tailored to model specific CME events, so the background solar wind solutions 
and resulting CME properties are not identical. Thus, the simulation and observational 
comparisons in this paper are, for the most part, qualitative. The charge state analyses in 
these simulations provide us with the full iron charge state distribution, so we can construct 
both the range of the distribution available in the in-situ data (Fe +6 to Fe +20 ) as well as 
the average iron charge states (Qp c ) f° r comparison to the spatial mapping in Figure [3J In 
Figure H] we plot the average iron charge state (Q-p e ) derived from the magnetic breakout 
model simulation (left panel; hereafter BM) and the flux cancellation model simulation 
(right panel; hereafter FC). The (Qf c ) contour maps are each made up of a 60-by-60 grid 
of Lagrangian fluid parcels where we have constructed {Qpe) = $Z (<7 n g/ Yl n q) fr° m the full 
distribution of iron charge states in each pixel. Representative flux contours are over-plotted 
showing the magnetic field cross sections of the flux rope ejecta, their sheath regions, and 
the surrounding solar wind structure. We note that the BM simulation results herein are 
presented with a flipped y-axis (compared to lLynch et al.ll201ll ) to ease the visual association 
with the spatial mapping results. 



3. Comparison of Observations and Simulation Results 

3.1. Large-Scale (Qfc) Structure 

The large-scale heliospheric structure of the average iron charge state signatures as- 
sociated with the May 21-23 successive ICMEs in Figure [3] have characteristics found in 
both the BM and FC MHD simulation results. MCI has only slightly enhanced average 
iron charge states of (Qfc) ~10.5 in the ACE measurements in the back half of the flux 
rope. MC2 has stronger enhancements, up to (Qfc) ~12, that are seen in both the STA 
and ACE trajectories through the flux rope. In both MCI and MC2 the location of the 
average iron enhancements are in the back half of the actual flux ropes, offset from their 
respective symmetry axes (flux rope center) , and also in the trailing solar wind flow following 
the rear boundaries of the flux rope field rotations. Looking at Figure HJ both the BM and 
FC simulations have a similar qualitative structure in terms of enhancements in the latter 
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half of the flux rope resulting from the contribution of flare-heated plasma being supplied 
by the eruptive flare reconnection jet. The magnitude of the charge state enhancements in 
the FC simulation are much larger, including (Qpe) values of >15 in the reconnection jet 
outflow and result in a shell of (Q-pe) ~ 12-13 enhancement encircling the flux rope's lower 
charge state core. In the BM simulation results there is also a clear enhancement preceding 
the flux rope in the topologically well defined sheath region due primarily to compressional 
heating of the density enhancements and in some part to the breakout reconnection heating 
during the eruption. This compares well with the charge state enhancement ahead of MC2 
which extends well beyond the flux rope boundaries. 

We can also look at Figure [3] in the context of the source region EUV dimmings of 
Figure[TJ For MCI, the charge state enhancements are minimal and follow the ICME flux 
rope. We see that in the MCI source region, the flare occurs on the eastern side of the 
eruption (dimming) which would have crossed 1AU beyond the STB spacecraft and so is 
not seen in the in-situ observations (i.e., the half of the MCI flux rope not included in the 
(Qfc) spatial map). For MC2 the bulk of the enhancements are in the western half of the 
event, which corresponds well to location of the flare in relation to the CME source region. 
The dimming associated with MC2 occurs rather strongly to the western side of the active 
region in which the flare occurs. The correspondence of the flare location and EUV dimming 
offset is consistent with the idea that the flare is indeed responsible for the CME heating 
that results in the enhanced charge states. 



3.2. Multispacecraft and Simulation Magnetic Field Structure 



Inspired by the general agreement of the large-scale average iron charge state signatures 
in the multispacecraft spatial mapping and the results from the MHD CME modeling, we 
now extend our comparison to a broader range of the iron charge state distributions sampled 
along synthetic spacecraft trajectories through the simulation results. The MCI and MC2 
flux ropes are each intersected by two of the three spacecraft and have the third spacecraft 
track either right at or just outside the flux rope boundary. The multispacecraft in-situ 
data yield six ICME intervals to be examined. In the top panel of Figure |5] (adapted from 
Mulligan et al .1120121 ) we plot the STEREO and ACE spacecraft trajectories through the MCI 
and MC2 flux rope ecliptic plane/ spatial map field rotations. T he MCI (MC2) ellipse colored 
yellow (blue) show the boundaries of the lMulligan et al.l ( 120121 ) multispacecraft elliptical flux 
rope model fits to the ICME flux rope interval. Figure [5] lower-left and lower-right panels 
show the BM and FC simulation snapshots of representative flux surfaces and the set of 
synthetic spacecraft trajectories that roughly correspond to the relative position and angle 
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of the STB, ACE, and STA tracks through MCI and MC2. For MCI (MC2) we chose the 
synthetic spacecraft trajectories at an angle of +22° (—12°) with respect to the simulation's 
circular flux rope cross section to qualitatively match the multispacecraft observing geometry. 
Thus, for each MHD simulation we define six model cuts to compare with the MCI and MC2 
data. For example, for the MCI interval we have the MCIstb, MCIace, and MCIsta in-situ 
data to be compared with the breakout model trajectories {BMg^ B , BM^ E , BMg^ A } and 
the flux cancellation trajectories {FCg^ B , FC^q E , FCg^ A }. 

Since the simulation trajectories represent spatial sampling through the BM and FC 
results and we are using the in-situ data time series, we want to ensure that the position 
of the simulation cut and time of the in-situ data correspond, as much as possible, to the 
same relative position with respect the magnetic flux ropes and the surrounding ICMEs. 
MCI and MC2 have di fferent flux rope orien tations (as seen in Figures [2] and of this paper 



and Figures 8 and 9 of Mulligan et al.ll2012l ). including axial field directions of opposite sign 
in the two events (seen largely as the sign of the spacecraft Bn component). Therefore, 
we have transformed the MHD simulation magnetic field data so that the sense of the field 
rotations and the direction of the axial component in the simulations match each of the 
MCI and MC2 events. In the axisymmetric MHD simulations, the (B r , Bg) field contain the 
azimuthal (twist) component of the flux rope ejecta while the B^ component represents the 
axial flux rope field. Multiplying the B^ component by —1 changes the orientation of the 
axial field but also changes the handedness of the flux rope unless the (B r , Bg) components 
are also multiplied by —1. 

Figures [6] and [7] plot the entire set of MCI and MC2 multispacecraft in-situ magnetic 
field and iron charge state data and the corresponding BM and FC simulation trajectories. 
The upper panels show the in-situ magnetic field components in each ICME interval, col- 
ored in the same manner as Figure HI i.e. B R (red), B T (green), B N (blue), B (black), and 
the fields from the MHD simulation trajectory sampling B r (red), Bg (green), B^ (blue), 
B (black). Obviously the simulation field magnitudes at ~12i? and the in-situ 1AU field 
magnitudes are very different. We have selected plot ranges to allow a quick visual com- 
parison to assess the simulation flux rope signatures with respect to the in-situ magnetic 
field. Here the spacecraft data are all plotted on the same ±20 nT range whereas the plot 
ranges for the fields along the BM simulation trajectories are ±0.02 G and along the FC 
simulation trajectories are ±0.08 G. Once again, and by construction, the agreement is quite 
reasonable. Both the BM and FC simulations capture the qualitative features of the flux 
rope rotation, at least in the two sampling trajectories that intersect the flux rope in each 
event: MCIstb, MCIace and MC2ace, MC2 ST a. This gives us confidence that the compar- 
ison of the simulation-derived ionic composition signatures and the in-situ ionic composition 
signatures is reasonable and can be used to both interpret the observations and evaluate the 
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simulations. 



3.3. Multispacecraft and Simulation Qp e Distributions 

The lower panels in Figures [6] and [7] plot the detailed iron charge state distributions of 
the in-situ data and from the MHD simulation trajectory sampling. For the Figure [6] MCI 
comparison, the ordering from top to bottom of STB, ACE, and STA represent cuts through 
the center, mid-range, and just missing the boundary of the magnetic flux rope. The in-situ 
Qpe distributions are consistently low throughout MCIstb, ranging primarily from +8 to 
+9. For MCIace the charge states have a decreasing profile starting with charge states of 
+9 to +12 and ending with charge states of mainly +8 and +9. Charge states in MCIsta 
range from +8 to +10 at the beginning and end of the event and drop slightly to +7 to +9 
at the center of the interval. 

Here, the BM simulation shows the general trend of the center of the flux rope having 
lower iron charge states: the peak of the distribution dips most strongly in BMg!^ B , more 
moderately in BM A( l E , and exhibit no variation in BMg!^ A . The background solar wind BM 
charge state distribution outside of the flux rope and sheath region is strongly peaked at 
+11, which is slightly higher than the no-enhancement levels seen during the ICME periods. 
The FC simulation charge states are more varied reflecting the intensity of the shell of flare 
heated material and corresponding charge state enhancement in the back of the flux rope 
structure. While the no-enhancement level of the FC distribution (+7 to +8) tends to be 
lower than the in-situ observations, for the simulation trajectories FCg^ B and FC A( l E we see 
a sharp transition to high iron charge states of +16 and a more gradual fall off of the level 
of enhancement in the trailing half. The ring structure in the enhancements is also visible in 
the FC cuts, with FC A( l E showing a narrow spike of +16 enhancement in the first half of the 
flux rope, but also as a weak signal in the FCg^ B and FC S ^ A plots. In the FCg!^ A trajectory 
there is also a remnant of the trailing half enhancement with the distribution peak shifting 
slightly higher to +8 to +9 and broadening to include trace amounts of +10 to +11. 

For the Figure [7] MC2 comparison, STB just misses the flux rope, while ACE and STA 
represent cuts through the center and mid-range impact angles. The in-situ Qp e distributions 
in MC2 have more structure and much higher levels of enhanced charge states. While 
MC2stb is still low (+8 to +9 peak), the distribution is broader with trace amounts of +11 
to +13. The latter half of the MC2 A ce distributions include a moderate amount of +11 to 
+13 and a small peak at +16. MC2sta exhibits the most elevated iron charge states with 
a bimodal structure including a strong peak at +15 to +16 and +8 to +9 at the center of 
the flux rope interval lasting ~8 hours. The bimodal iron charge state structure is common 
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feature of ICME ev ents with flare heated composition enhancements (IGruesbeck et al.ll2011 



Gilbert et al.ll2012j ). 



The BM simulation charge state data for the MC2 comparisons show the same trend as 

(2) 

the MCI cuts in terms of distance from the flux rope center; BM A( ^ E shows the largest central 

(2) (2) 

lower charge state dip, BM S1 ; A shows a shallower dip, and BM gi ; B shows very little variation. 
Again, the baseline BM distribution peaks are slightly higher than the non-enhanced periods 
of the MC2 distributions. The iron charge state distributions in the FC simulation match 
features of the MC2 in-situ data reasonably well for the two flux rope intersections. Here 

(2) 

the MC2ace enhancements and the FC A £ E distributions appear in the same trailing region 
of the flux ropes and the gradual return of the iron charge state distribution from +16 
to lower enhancement levels are similar. The FClv^ sampling trajectory stays in the flare 
heated enhancement ring longer and therefore shows the most extended hot iron +16 peak 
in the center of the magnetic field rotation signatures. The relative duration of this central 

(2) 

enhancement is qualitatively similar to the MC2sta data. The FCgrf A enhancement region 
also has a bimodal distribution, with the low charge state peak at +8 in the first half of the 
region then elevating to higher charge states about one third of the way through the event 
and staying elevated until about two thirds of the way through the event. This structure is 
qualitatively similar to what is seen in the STA observations. 



4. Discussion 



In our comparisons of Qp e derived from MHD simulations and multispacecraft ICME 
observations we find several similarities and differences. In Figures E] and [7J the cuts are 
chosen to approximate the magnetic field rotations seen at each spacecraft. The ambient 
level of Qyg in the data is lower than the BM and higher than the FC model. The specific 
charge state enhancements for MCI are not well matched to the observations. For example, 
the STB observations of MCI, which crossed very close to the center of the flux rope, show 
Q-pe ranging primarily from values of +7 to +12. The BM results for a similar cut have 
charge states ranging from +8 to +14 with a distinctive dip at the flux rope center, while 
the FC results have values that begin at +7 to +8 which then jump to +16 at the location 
of the flare reconnection heated material. There are, however, some similarities between 
MC2 and the FC model, particularly for the ACE and STA cuts. In the ACE cut, the data 
show an abrupt jump to charge states of +16 about halfway through the flux rope, which is 
similar to the jump seen in the FC cut and extends to the end of the flux rope. At STA the 
charge state increase occurs about 1/3 of the way through the event and is mostly confined 
to the center third of the event, in both the data and the FC model results. 
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In the global view of the spacecraft data (Figure [3]) we find that the bulk of the charge 
state enhancements are seen in the back half and surrounding the flux ropes. For MCI 
the enhancements are relatively low and appear in the back half of the flux rope and the 
region following the flux rop e. For MC2, the charge st ate levels (Qp c ) > +12 are consistent 
with ICME related material (ILepri fc Zurbuchenl 12004] ) . The center of the flux rope contains 
charge states that are similar to ambient solar wind values, while the enhancements seem to 
surround the flux rope, being located primarily in the trailing portion of the flux rope, along 
the western (STA) flank, and in the sheath region ahead of the flux rope structure. In the 
two simulations (Figure H]) we find that the center of the flux ropes are "colder" than the 
ambient solar wind values with charge state enhancements present surrounding and trailing 
this colder region. Interestingly, in the BM we see enhancements preceding the flux rope 
that are very similar to the enhancements seen the sheath region of MC2. These two flux 
ropes, particularly MC2, bear an encouraging qualitative resemblance to the model results. 

In both events and the two MHD models the charge state enhancements are present 
most strongly in the second half of the flux rope and in the trailing material. This result 
is surprising as previous studies have found that charge state data is more likely to be 
enhan ced near the center of the event, and in particular to be cospatial with the flux rope 



regionf lHenke et al.lll998l bOOll iLvnch et al.ll2003l : iRichardson fc Canel 12004] ; iReinardl 12005 



20081 ). One possible explanation is that these events illustrate a different or less efficient 
mechanism that is less common and t hus do not stand out in the statistical sampling of 
Henke et al.l ( 120011 ) and IReinardl ( 120081 ). Since the associated fla res were only small B-class 
flares and the enhanced charge states only slightly exceed the ILepri &: Zurbuchenl (120041 ) 
threshold of +12 required for enhanced (Qf c ), it is possible that in these events the small 
associated flares did not have sufficient energy to heat the center of the flux rope. An event 
associated with a larger M or X-class flare, on the other hand, could experience a more 
explosive eruption in which the reconnection starts at the center of the event and moves 
outwards causing a larger flare and a CME with the strongest heating in the center. In that 
case, the much larger enhanced charge states in the center of the event would overwhelm the 
post-eruption flare arcade related heating that we may be seeing here. 

On the other hand, the previous statistical studies may have averaged over the fine 
structure of the ionic composition enhancements. These studies were based on event averages, 
which would blur the boundary between enhancements present in an internal ring around 
the flux rope center and the less enhanced central region, such as that seen in these events 
(Figure [3]) and the charge states derived from MHD models (Figure H]). In addition, the 
previous studies were not designed to examine asymmetries in the structure, so enhancement 
offsets, such as seen in Figure E] would be averaged out. Additional multispacecraft studies, 
particularly of CMEs with a large associated flares, would help resolve this issue. In the 
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meantime, more precise statistical studies are planned, focusing on spatial position within 
the ICME and relative location of the flare/AR site in order to shed some light on this 
question. 



Recent numerical MHD modeling of sympathetic eruptions by iTorok et al.l (120111 ) pro- 
vide a possible scenario to explain the the complex spatial structure of the interplanetary 
ionic co mposition enhancem ents observed in these events. Here we refer the reader to Fig- 



ure 3 of ITorok et al.l ( 120111 ) and specifically the location of the current sheets and resulting 
reconnection regions associated with the sympathetic eruptions of their second and third flux 
ropes. The authors simulate the consecutive, linked eruptions of three flux ropes to mimic 
the August 1, 2010 series of filament eruptions. The consecutive May 19-20, 2007 eruptions 
resulti ng in our in-s i tu obs ervations of the MCI and MC2 ICMEs bear a striking resemblance 



to the ITorok et al.l ( 120111 ) scenario's second and third filament /flux rope eruptions. While 
the standard eruptive flare (Forbes, 2000 and references therein) current sheet underneath 
each flux rope eruption supplies the bulk flare heating and enhanced Q? e signatures to the 
latter portion of the flux rope ejecta, in the sympathetic eruption scenario, their second 
eruption current sheet (May 19 CME — > MCI) also acts as the overlying magnetic breakout 
reconnection for the third flux rope eruption (May 20 CME — > MC2). Here the May 19 
CME flare-heated QFe enhancements occur continuously right up until the eruption of the 
May 20 CME, providing a natural explanation for the multispacecraft in-situ charge states 
signatures that seems to link MCI and MC2 in our Figure [3] spatial mapping results: both 
a moderate flare-heated Qp e enhancement in the trailing portion of the flux rope interiors 
and an extended, continuous enhancement in the solar wind period between the two flux 
ropes. Taking into account the ionic charge state signatures actually allows us to unravel 
and understand the complex field and plasma signatures in the various spacecraft data sets 
during our "dual flux rope" composite/compound event. 

It is interesting that the lower half of MC2 (STB measurements) do not show a similar, 
symmetrical enhancement. As explained earlier, however, we see that the flare is not sym- 
metrically located under the CME (unlike our simulations). If we assume that the dimming 
region represents the eruption site, the active region and flare occur to the west of the MC2 
source region. This is consistent with heating of the west side of the CME as seen in in-situ 
data. In the MCI eruption, on the other hand, the active region and flare occur to the east 
of the source region dimming and it is thus possible that the eastern half of the MC, which 
is unseen, has more enhancements. 

Finally, we note that ionic and elemental composition can be measured in the corona 
"directly" via techniques utilizing UV spectroscopy. The SOHO/UVCS observations, typi- 
cally combined with other multi-wavelength emission data, can be used to infer the coronal 
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densities and temperatures derived from various line ratios as well as estimating the amount 
of fl are and CME-related heating during the early phases of the CME eruption 5,2-R^ (see, 



Akmal et all 120011 iKo et all 120031 : Bemporad et all 120071 iLee et all 1200a lLandi et al 



20101 ). While these analyses are beyond the scope of the current paper, we acknowledge that 



bridging the gap between coronal spectroscopy and in-situ composition measurements, both 
in the solar wind and in CMEs, is an important arena for future research. 



5. Conclusions 

We have described an effort to understand the full global structure of ICMEs from the 
Sun to the Ea rth by combining a novel method to derive ionic charge state distributions from 
MHD models ( Lynch et alJboil ) with newly developed, groundbreaking multi spacecraft dat 



analy sis techniques that allow us to construct the global structure of ICMEs ([Mulligan et al. 



20121 ). We compare the iron charge states derived from idealized MHD simulations with 



observations of the May 21-23, 2007 ICMEs. 

The large scale comparisons (Figures 3 and 4) show similarities between the Qp e derived 
from MHD simulation results and the spatial mapping results derived from STEREO and 
ACE observations of the May 2007 events. In particular, we find that the charge state en- 
hancements tend to occur at the back end and surrounding the flux rope material. Both the 
BM and the FC models have similar behavior with enhanced charge states surrounding and 
following the colder flux rope center. Enhancements seen beyond the flux rope in the data, 
particularly ahead of the MC2 ICME, are well matched to the BM results, suggesting MCI 
eruption's eruptive flare reconnection was simultaneously acting as breakout-like reconnec- 
tion during the sympathetic eruption scenario resulting from the source region's multipolar 
topology. 

More detailed time series comparisons between the spacecraft observations and cuts 
through the MHD model results (Figures 6 and 7) show that the ambient solar wind charge 
state levels are lower than predicted by the BM and higher than predicted by the FC model 
simulations. Quantitative agreement between the observed flux rope iron charge state dis- 
tributions are not well captured in the modeling results, particularly for MCI. However, 
the FC model does remarkably well at qualitatively matching the increased charge states 
observed in MC2, including both the Qpe distribution broadening in the ACE data and the 
bimodal Qf g signature in the STA data. Therefore, we find aspects of both the BM and 
FC simulation- derived charges states present in the 2007 May 21-23 multispacecraft ICME 
observations. 
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Overall, our analysis and results provide the most comprehensive view to date of the 
global structure of enhanced composition within ICMEs. More investigation is needed to 
determine whether the structures observed in this case study are common; in particular, the 
appearance of enhanced charge states surrounding the flux ropes is at odds with previous 
statistical results indicating that enhanced charge states are more common at the flux rope 
center. Finally, we show that the ability to derive charge state information from MHD 
models provides an important and potentially very useful constraint for these models. 

A.A.R. and T.M. acknowledge support from NASA SR&T NNX08AH54G. B.J.L. ac- 
knowledges support of AFOSR YIP FA9550-11-1-0048 and NASA HTP NNX11AJ65G. Sup- 
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Fig. 1. — EUVI 171 A base (top) and difference (bottom) images showing changes in the 
source region following CME eruptions on May 19 (left) and May 20 (right). These dimmings 
indicate where the majority of the CME plasma originated. In the case of the May 19 (MCI) 
event, the eruption took place primarily to the northwest of the active region, while in the 
May 20 (MC2) event the eruption was located to the southeast of the active region. 
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Fig. 2.— STEREO-B (left), ACE (middle), and STEREO-A (right) spacecraft data for 
2007 May 21-24 (DOY 141-144). From top to bottom shows v p , n p , T p , and the empirical 
"expected" temperature T exp (green), distribution of iron charge states (from +6 to +20), 
and the RTN magnetic field components B R (red), B T (green), £? N (blue), and B magnitude 
(black). The event interval MCI (MC2) is indicated as vertical yellow (light blue) lines with 
trajectories that intersect the ICME flux rope shown as the shaded yellow (light blue) bars. 
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Fig. 3. — Application of the Mulligan et al.l ( 12012I ) spatial mapping interpolation technique 
to the average iron charge state (Qp e ) data for the 2007 May 21-24 interval. The gray lines 
show, from top to bottom, the STEREO-A, ACE, and STEREO-B spacecraft tracks, i.e. the 
position of the each spacecraft's measured data. The vector field arrows show the magnetic 
field in the ecliptic plane during the MCI and MC2 regions indicated in Figure [2j 
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Fig. 4. — Snapshot of the average iron charge state (Qpe) spatial distribution derived from 
2D MHD simulations of CME initiation: Left panel, the breakout model (BM); Right panel, 
the flux cancellation (FC) model. The simulation flux ropes are at ~12i? in their respective 
domains but the heavy ion charge states have already frozen-in. Snapshots were selected so 
that each flux rope was at or past 12 R & . See iLynch et al.l (120111 ) for further details. 
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Fig. 5. — Top panel plots a portion of the multispacecraft spatial m apping tracks and th e 
ecliptic plane field rotations for the two ICME intervals (adapted from Mulligan et al.ll2012l ). 
Lower left panel shows the synthetic spacecraft trajectory cuts through the breakout model 
simulation corresponding to MCI (BM^) and MC2 (BM*- 2 * 1 ). Lower right panel shows the 
trajectories through the flux cancellation simulation (FC*- 1 - 1 , FC^ 2 -*). 
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Fig. 6. — Comparison between the MCI multispacecraft data and the MCI simulation 
trajectory samples shown in Figure The top, middle, bottom plots show the intervals for 
STEREO-B, ACE, and STEREO- A respectively. Each plot shows the observed magnetic 
field (upper panel), the extended QF e distribution (lower panel), and the corresponding 
quantities in the breakout model (BM) and flux cancellation (FC) simulations. In order to 
facilitate the comparison between the observations and simulation data we have transformed 
the simulation magnetic field profiles to match the sign of the axial flux rope field and sense 
of rotation; see text for details. 
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Fig. 7. — Comparison between the MC2 multispacecraft data and the MC2 simulation 
trajectory samples shown in Figure [5] in the same format as Figure [6j 



